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Small x physics and hard QCD processes at LHC 



1. Itroduction 



The so-called small x regime of QCD is the kinematic region, where the characteristic hard 
scale of the process \i 2 ~ p 2 ~ Mj = M 2 + p\, ( M and p\ are the mass and the transverse 
momentum of the final state) is large as compared to the Aqcd but ji. is much less than the total 
c.m.s. energy y/s of the process: Aqcd <C /J. <C v^- 

In this sense, the HERA was the first small x machine, and the LHC is more of a small x 
collider. Typical x values probed at the LHC in the central rapidity region are almost two orders 
of magnitude smaller than x values probed at the HERA at the same scale. Hence, the small x 
corrections start being relevant even for a final state with a characteristic electroweak scale M ~ 100 
GeV. 

It means the pQCD expansion any observable quantity in a s contains large coefficients 
(In" (S/M 2 )) ~ (ln"(l/x)) (besides the usual renorm group ones (ln"(/i 2 /Ag CZ) )). The resumma- 
tion of these terms (a 4 (ln(l/x))" ~ 1 at x — > 0) in the framework of the Balitsky-Fadin-Kuraev- 
Lipatov (BFKL) theory §\§ results in the so called unintegrated gluon distribution ^(x,kj). The 
unintegrated gluon distribution determines the probability to find a gluon carrying the longitu- 
dinal momentum fraction x and transverse momentum kj. This generalized factorization is called 
"k T — factorization" Q|. If the terms proportional to a"ln n (jx 2 / A 2 QCD ) and a s ln n {ji 2 / A 2 Q CD )ln n (\ / x) 
are also resummed, then the unintegrated gluon distribution function depends also on the probing 
scale fl. This quantity depends on more degrees of freedom than the usual collinear partem density, 
and is therefore less constrained by the experimental data. Various approaches to model the unin- 
tegrated gluon distribution have been proposed. One such approach, valid for both small and large 
x, has been developed by Ciafaloni, Catani, Fiorani and Marchesini, and is known as the CCFM 
model It introduces angular ordering of emissions to correctly treat gluon coherence effects. 
In the limit of asymptotic energies, it is almost equivalent to BFKL but also similar to the 
collinear (DGLAP) evolution for large x and high \i 2 . The resulting unintegrated gluon distribution 
functions depend on two scales, the additional scale q being a variable related to the maximum 
angle allowed in the emission. 

The BFKL evolution equation predicts rapid growth of gluon density (~ x~ A , where 1 + A is 
the intercept of so-called hard BFKL Pomeron). However it is clear that this growth cannot continue 
for ever, because it would violate the unitarity constraint [g]. Consequently, the parton evolution 
dynamics must change at some point, and new phenomenon must come into play. Indeed as the 
gluon density increases, non-linear parton interactions are expected to become more and more 
important, resulting eventually in the slowdown of the parton density growth (known as "saturation 
effect") ^ ||]. The underlying physics can be described by the non-linear Balitsky-Kovchegov 
(BK) equation [g]. These nonlinear interactions lead to an equilibrium-like system of partons with 
some definite value of the average transverse momentum kj and the corresponding saturation scale 
Q s (x). This equilibrium-like system is the so called Color Glass Condesate (CGC) [{J. Since the 
saturation scale increases with decreasing of x: Q 2 (x,A) ~ x~ x A 8 (A is an atomic number ) with 
A ~ 0.3, 8 ~ 1/3 one may expect that the saturation effect will be more clear at LHC energies. 
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2. Unintegrated parton distributions (uPDF or TMD) 

The basic dynamical quantity in the small x physics is transverse-momentum-dependent (TMD) 
(kr -dependent) or unintegrated parton distribution (uPDF) £/(x,kj,iJ. 2 ). For example to calculate 
the cross sections of photoproduction process the uPDF si (x,k|,/i 2 ) has to be convoluted with the 
relevant partonic cross section a yg : 

a = JjJ dk 2 T o rg (x/z,k 2 T ,ii 2 )^(x,k 2 T ,ii 2 ). (2.1) 

For the uPDF there is no unique definition, and as a cosequence it is for the phenomenology of these 
quantities very important to identify uPDF which are used in description of high energy processes. 
For a general introduction to small x physics and the small x evolution equations, as well as tools 
for calculation in terms of MC programs, we refer to the reviews 

During roughly the last decade, there has been steady progress toward a better understanding 
of the kj -factorization (high energy factorization) and the uPDF (for example [[K]]. Workshop on 
Transvere Momentum Distributions (TMD 2010), which was held in Trento (Italy), was dedicated 



to the recent developments in small x physics, based on the kj -factorization and the uPDF [11]. 



Recently the definition for the TMDs determined by the requirement of factorization, maxi- 



mal universality and internal consistency have been done by Collins Q12|]. The results obtained in 



previous works are reduced to the following: &r(TMD)-factorization is valid in 

• Back-to-back hadron or jet production in e + e~ -annihilation, 

• Drell-Yan process (P a +Pb -> (f,W/Z) +X), 

• Semi-inclusive DIS {e + P ^ e + h + X). 

In hadroproduction of back-to-back jets or hadrons ( h\ + h% — I H\ +H2 +X) TMD-factorization 
is problematic. 

For expample, partonic picture gives the following qj -dependent hadronic tensor for DY cross 
section: 

W^=L f \H f (Q;iiR)r (2.2) 
J d 2 kiTd 2 k2T^f/p l (xi,kiT^R^i)£/j / p 2 (x2,k2T^R^2)8(kiT + k2T-qT)+Y(Q,q T )^ 

The hard part Hf(Q; jj, R ) is calculable to arbitrary order in a s , /j,r - the renormalization scale. 
The term Y(Q,qr) describes the matching to large qj, where the approximations of TMD-factorization 
break down. The scales £j , £ 2 are related to the regulation of light-cone divergences and £j §2 = Q 4 . 
The soft factors connected with soft gluons are contained in the definitions of the TMDs, which 
cannot be predicted from the theory and must be fitted to data. 

3. The kj -factorization approach in hadroproduction 

The kj -factorization approach in hadroproduction is based on the work by Catani, Ciafaloni 
and Hautman (CCH) (see [[J) The factorization formula for /7/7-collision in physical gauge (nA = 
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0,n^ =aPf + Wf )is 

° = W 2 / d2klT J~x^~J d2k2T /^^(■ x i' k ir)^(P/( x i- x 2),kir,k 2 r)^(^2,k2r), (3.1) 

where p = AM 2 /s, M is the invariant mass of heavy quark, and & are the unintegrated gluon 
distributions, definded by the BFKL equation: 

^(x,k;Q 2 ) = U(l-x)S(k 2 -Q 2 ) + (3.2) 

+a s f^f ^(x/z^ + q-QD-eik-q^ix/z^Ql)}, 
J nq z J z 

were a s = a, s N c /%. It means that the rapidity divergencies are cut off since there are an implicit 
cuts in the BFKL formalism. Effectively one introduces a cuts £1 , £ 2 , and then sets £j = x\ , £ 2 = x 2 
in (3.1). 

The declaration in CCH [0] that is defined via the BFKL equation (3.2) means that the BFKL 



unintegrated gluon distribution reduces to the dipole gluon distribution [13]. The connections 



between different uPDF recently were analysed in [14]. 

The procedure for resumming inclusive hard cross-sections at the leading non-trivial order 
through kj -factorization was used for an increasing number of processes: photoproduction ones, 
DIS ones, DY and vector boson production, direct photon production, gluonic Higgs production 



both in the point-like limit, and for finite top mass m t . Please look, for example [15] 



The hadroproduction of heavy quarks was considered in [[If]] and recently in []17p. In last paper 
it was shown that when the coupling runs the dramatic enhancements seen at fixed coupling, due 
to infrared singularities in the partonic cross sections, are substantially reduced, to the extent that 
they are largely accounted for by the usual NLO and NNLO perturbative corrections. It was found 
that resummation modifies the B production cross section, at the LHC by at most 15%, but that the 
enhancement of gluonic W production may be as large 50% at large rapidities. 

In our previous papers we have used the kj -factorization approach to describe experimental 
data on: 

• heavy quark photo- and electroproduction at HERA 

• J/y production in photo- and electroproduction at HERA with taking into account the color 
singlet and color octet states. 

• D*,D* + jet, D* + 2 jet photoproduction and D* production in DIS 

• charm contribution to the structure function F 2 C '(x, Q 2 ),F[,Fl 

• B-meson and bb pair production at the Tevatron 

• charm, beauty, D* and J j\\f production in two-photon collisions at LEP2 

• Higgs production at the Tevatron and LHC 

• prompt photon production at the HERA and Tevatron 

• W/Z production at the Tevatron 



4 



Small x physics and hard QCD processes at LHC 



Here I want to present the results of ft-quark and J / production at the LHC [18, |l^]in com- 
parison with first experimental data obtained by ATLAS, CMS and LHCb Collaborations. The 



description of prompt photon production and DY lepton pairs was done by M. Malyshev [|20|]. 



4. Ingredients of our fc r -factorization numerical calculations 



To calculate the cross section of any physical process in the kj— factorization approach accord- 
ing to the formula (3.1) the partonic cross section a has to be taken off mass shell (kr -dependent) 
and the polarization density matrix of initial gluons has to be taken in the so called BFKL form 1 : 



'e"e* v 



(4.1) 



Concerning the uPDF in a proton, we used two different sets. First of them is the KMR 
one. The KMR approach represent an approximate treatment of the parton evolution mainly based 
on the DGLAP equation and incorpotating the BFKL effects at the last step of the parton ladder 
only, in the form of the properly defined Sudakov formfactors T q (kj,jj. 2 ) and 7^(ky, jU ), including 



logarithmic loop corrections [21] 



x dz 



^(x ! k 2 r , M 2 ) = r ? (k|, M 2 )^x 



(4.2) 



l 

X J ' dZ 

X 



(x,k 2 ^ 2 ) = T g (k 2 ,^ 



aJk 



2n 



■z) 



(4.3) 



where ©-functions imply the angular-ordering constraint A = n/(n+kr) specifically to the last 
evolution step (to regulate the soft gluon singularities). For other evolution steps the strong or- 
dering in transverse momentum within DGLAP equation automatically ensures angular ordering. 
T a (kj,iJ. 2 ) - the probability of evolving from k 2 to pL 2 without parton emission. r o (kj-,/i 2 ) = 1 at 
ky > pL 2 . Such definition of the srf a {x,k\ ,pL 2 ) is correct for k 2 > /I 2 only, where /io ~ 1 GeV is 
the minimum scale for which DGLAP evolution of the collinear parton densities is valid. 



We use the last version of KMRW uPDF obtained from DGLAP equations [|22[]. In this case 
(a(x,n 2 ) = xG or a(x,/x 2 ) = xq) the normalization condition 



a(x,n 2 ) = J £/ a (x,kj,n 2 )dkY 



(4.4) 



The problem of choicing of proper gauge will be discussed in more details in Sec. 5. 
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Figure 1: Feynman diagrams for the pp — > QQX process. 



is satisfied, if 



2 ,,2-1 



< (x, k T , | k 2 <M 2 = a (x, Mo ) T a (At , H 



(4.5) 



where T a (fj^,fi 2 ) are the quark and gluon Sudakov form factors. Then the uPDF £/ a (x,kj,iJ. 2 ) is 
defined in all kf. region. 

Another uPDF was obtained using the CCFM evolution equation. The CCFM evolution equa- 
tion has been solved numerically using a Monte-Carlo method [ |23| ] 

According to the CCFM evolution equation the emission of gluons during the initial cascade is 
only allowed in an angular-ordered region of phase space. The maximum allowed angle E related 
to the hard quark box sets the scale /i: p} = S + Qj (= jxj). 

The unintegrated gluon distribution are determined by a convolution of the non-perturbative start- 
ing distribution s/q{x) and CCFM evolution denoted by £/(x,kj,n 2 ): 



X£/(x,kj,jJ,' 1 



z z 



(4.6) 



where 



xM>(x) = Nx P0 (\-x) P[ exp(-\4/k%). 
The parameters were determined in the fit to F2 data. 



(4.7) 



5. Heavy quark production in pp-interaction 



The hard partonic subprocess g*g* — > QQ is described by the Feynman's diagrams presented 
in Fig. 1. We used Sudakov decopmosition for the momenta of heavy quarks and the initial gluons: 



Pi = aiPi+piP 2 +Pi±, h = aPi+k l± , k 2 =pP 2 + k 2 ±,p 2 l 
in the center of mass frame of colliding particles 



Pi 



■ k 2 T , where 
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Pi = (E,0,0,E), P 2 = (E,0,0,-E), E = ^s/2, Pf = P% = 0, (P l P 2 )=s/2. 
Sudakov variables are 

M \T , *x M 2 T , *x „ M\ T , _ „ M 2 T , * N 

«i = — —expiyi), a 2 = — — exp(j 2 ), ft = — exp^-yj), j8 2 = — ^exp(-y 2 ), 
V s v 5, V 5, 

ki T +k 2T = pit + p 2 t, a = a l + a 2 , /} = j8i+j8 2 - 

To guarantee gauge invariance, the process with off-shell incoming particles has to be embed- 
ded into the scaterring of on-shell particles. The second row of Fig. 1 includes non-factorizing 
diagrams which are factorized only in the sum. To make this factorization one can sum up these 
diagrams with the last diagram in the first row leading to one diagram with an effective Lipatov 
vertex by working in Feynman gauge [ 24 ] : 

r (*„te) - ™ - - (hr - hrV) , (5.1) 

s \ ps as J 

where t\=k\ = -kj T , t 2 = kj = -k 2r , Mj = s+ (k\ T + k 2r ). 

This vertex obeys the Ward identity: r^(k\,k 2 )ki^ = 0. Then the last five diagrams in Fig. 1 are 
replaced by one diagram (Fig. 2). By neglecting the exchanged momentum in the coupling of 
gluons to incoming particles, we get an eikonal vertex which does not depend on the spin of the 
particle: 

u(X[ t PL-k 1 )'fu(X 1 ,P l ) ^2P»S KM . (5.2) 

Then it is possible to remove the external particle lines and attach so-called "non-sense" polariza- 
tion to the incoming gluons: 



; = V2P?/V~s, < = V2P^/V~s. (5.3) 



Instead of Feynman gauge, one can choose an appropriate axial gauge n • A = (n? = aPf + bP%). 
The contraction of the eikonal coupling with the gluon polarization in this gauge 

41 (*) = S, V + + - (5.4) 
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then reads 

Pf 4v (*i ) = kirv/cc, P$d$ (k 2 ) = k 2Tv /p . (5.5) 

In such a physical gauge the non-factortizing diagrammes vanish since the direct connection of two 
eikonal couplings gives P^d^yP^ = 0. It means the Lipatov vertex is to be replaced by the usual 
three gluon vertex. Then we can use the following matrix elements according to the diagrams in 
Fig. 1: 

Mi = M(^ 1 )(-/gy^)£ M (fc 1 )/ ^ 1 _^ 2 + _^ 2 (-/gy v )£ v fe)v( j p 2 ), (5.6) 
M 2 = u( Pl )(-igf)e v (k 2 )i ^Zllf M M 2 Hg^) g M (*i HP2), (5-7) 
M 3 = u{ Pl )C^ x {-k u -k 2 M +k 2 ) 8 \ {k ^ 2 2) Yxv(P2), (5.8) 

where 

C^ x {k u k 2 M)=i((k2-h) x g^ + {h-k 2 fg vX + {k l -hYg^) (5.9) 
is the standard three gluon vertex. 



6. Numerical results 



Recently we have demonstrated reasonable agreement between the kj -factorization predic- 



tions and the Tevatron data on the ^-quarks, bb di-jets, B + - and D-mesons Q25Q. Based on these 



results, here we give here analysis of the CMS [26, 27, 28] and LHCb [29] data in the framework 



of the kj -factorization approach. We produced the relevant numerical calculations in two ways: 

• We performed analytical parton-level calculations (which are labeled as LZ). 

• The measured cross sections of heavy quark production was compared also with the predic- 
tions of full hadron level Monte Carlo event generator CASCADE. 

In our numerical calculations we have used three different sets, namely the CCFM AO (BO) 
and the KMR ones. The difference between AO and BO sets is connected with the different val- 
ues of soft cut and width of the intrinsic distribution. A reasonable description of the F 2 
data can be achieved by both these sets. For the input, we have used the standard MSTW'2008 
(LO) [@ (in LZ calculations) and the MRST 99 [||] (in CASCADE) sets. The unintegrated 



gluon distributions depend on the renormalization and factorization scales /i« and jif. We set 
Mi? = m Q + (Pir + Kr)/^' = ^ + Qt> wnere Qt is the transverse momentum of the initial 
off-shell gluon pair, m c = 1.4 ±0.1 GeV, m h = 4.75 ±0.25 GeV. We use the LO formula for the 
coupling a s (jUjj) with «/ = 4 active quark flavors at Aqcd = 200 MeV, such that a 4 (M|) = 0.1232. 



We begin the discussion by presenting our results for the muons originating from the semilep- 
tonic decays of the b quarks. The CMS collaboration has measured the transverse momentum and 
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Figure 3: The pseudo- rapidity distributions ofmuons arising from the semileptonic decays of beauty quarks. 
The first column shows the LZ numerical results while the second one depicts the CASCADE predictions. 
The solid, dashed and dash-dotted, dotted histograms correspond to the results obtained with the CCFM AO, 
BO and KMR unintegrated gluon densities. The experimental data are from CMS 

the pseudorapidity distributions of muons from ^-decays. The measurements have been performed 
in the kinematic range pj>6 GeV and \ r\^\ < 2.1 at the total center-of-mass energy y/s = 7 TeV. 
To produce muons from ft-quarks, we first convert ^-quarks into B-mesons using the Peterson 
fragmentation function with default value Ef, = 0.006 and then simulate their semileptonic decay 
according to the standard electroweak theory taking into account the decays b — > pL as well as the 
cascade decay b — > c — > pL. 

The results of our calcilations are shown in Figs. 3 - 8 in the comparison with the LHC 
data (see fli8| ] for more details). We obtain a good description of the data when using the CCFM- 
evolved (namely, AO) gluon distribution in LZ calculations. The shape and absolute normalization 
of measured ^-flavored hadron cross sections at forward rapidities are reproduced well (see Fig. 6). 
The KMR and CCFM B0 predictions are somewhat below the data. In contrast with b hadron and 
decay muon cross sections, the results for inclusive Z?-jet production based on the CCFM and KMR 
gluons are very similar to each other anda reasonable description of the data is obtained by all 
unintegrated gluon distributions under consideration. 

The CASCADE predictions tend to lie slightly below the LZ ones and are rather close to the 
MC@NLO calculations (not shown). The observed difference between the LZ and CASCADE 
connects with the missing partem shower effects in the LZ evaluations. We have checked addition- 
aly that the LZ and CASCADE predictions coincide at parton level.. 

Figs. 5 and 8 show the role of fragmentation and off-shelness effects in our calculations. 

7. Quarkonium production in the kj— factorization approach 

The production of prompt J / i/A(T)-mesons in /^-collisions can proceed via either direct gluon- 
gluon fusion or the production of P-wave states XciXb) and 5-wave state y' followed by their 
radiative decays Xc(Xb) ■VvOO + 7- In the kj— factorization approach the direct mechanism 
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Figure 4: The transverse momentum distributions ofmuons arising from the semileptonic decays of beauty 
quarks. The first column shows the LZ numerical results while the second one depicts the CASCADE predic- 
tions. Notation of all histograms is the same as on previous slide. The experimental data are from CMS /p6|/. 




Figure 5: The dependence of our predictions on the fragmentation scheme. The solid, dashed and dash- 
dotted histograms correspond to the results obtained using the Peterson fragmentation function with £y = 
0.006, £y = 0.003 and the non-perturbative fragmentation functions respectively. We use CCFM (AO) gluon 
density for illustration. The experimental data are from CMS /^/. 



corresponds to the partonic subprocess g* +g* — > J / yQO +g. The production of P-wave mesons 
is given by g* + g* — > XciXb), an d there is no emittion of any additional gluons. The feed-down 
contribution from 5-wave state y' is described by g* + g* — > y' + g. 

The cross sections charmonium states depend on the renormalization and factorization scales }Xr 
and /if. We set /i^ = m 2 + pf- and jtij; = s + Qj, where Qj is the transverse momentum of ini- 
tial off-shell gluon pair. Following to PDG [32], we set = 3.097 GeV, m xc \ =3.511 GeV, 
m XC 2 = 3.556 GeV, m^i = 3.686 GeV and use the LO formula for the coupling constant a s {jl 2 ) 
with nf = 4 quark flavours at Aq C d = 200 Mev, such that a(M|) = 0.1232. 

The charmonium wave functions are taken to be equal to \Rj / w (0)\ 2 / An = 0.0876 GeV 3 , \R' x {0)\ 2 = 
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Figure 6: The pseudo rapidity distributions of b-flavored hadrons at LHC. LZ results with the CCFM AO, 
BO and KMR uPDF. The experimental data are from LHCb /^/. 
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Figure 7: The differential cross sections do j dydpj of inclusive b-jet production integrated over the speci- 
fied y intervals. The experimental data are from CMS / 



0.075 GeV 5 , \R ¥ >(0)\ 2 /4n = 0.0391 GeV 3 and the following branching fractions are used B(Xd -> 
J/Y+ y) = 0.356, B(^ c2 ^ J/W+7) = 0.202, = 0.561 and B(J/y^ jU + M~) = 
0.0593. Since the branching fraction for XcO ~^ J /W+ /decay is more than an order of magnitide 
smaller than for % c \ and % c %, we neglect its contribution to J/\j/ production. As y' — > J/y + X 
decay matrix elements are unknown, these events were generated according to the phase space. 



Comparison the results of our calculations with the CMS [33], ATLAS [34] and LHCb [35] 



data are shown in Figs. 9 - 11 [19]. We see that the taking into account sole direct production 



11 



Small x physics and hard QCD processes at LHC 




Figure 8: Importance of non-zero kj of incoming gluons. Dotted histograms - the results obtained without 
the virtualities gluons and with k|, < jx^ in matrix element. The CMS data /|2^/. 

is not sufficient to descibe the LHC data. We have obtained a good overall agreement between our 
predictions and the data when summing up the direct and feed-down contributions. The dependence 
of our numerical results on the uPDF is rather weak and the CCFM and KMR predictions are 
practically coincide. The difference between them can be observed at small pj or at large rapidities 
probed at the LHCb measerements. We have evaluated the polarizations parameters of prompt J / y 
mesons in the kinematical region of CMS, ATLAS and LHCb measurements in the Collins-Soper 
and helicity frames (see [|l9|]). We have took into account the contributions from the direcct and 
feed-dowm mechanisms. The qualitative predictions for the J/y meson polarization are stable 
with respect to variations in the model parameters. Therefore future price measurements of the 
polarizations parameters of the J/y mesons at the LHC will play crucial role in discriminating the 
different theoretical approaches. 

8. Conclusions 

In the present time there is steady progress toward a better understanding of the kj -factorization 
(high energy factorization) and the uPDF (TMD). 

We have described the first exp. data of Z?-quark and J/y production at LHC in the kj -factorization 
approach. We have obtained reasonable agreement of our calculations and the first experimental 
data taken by the CMS and ATLAS Collaborations. 

The dependence of our predictions on the uPDF appears at small transverse momenta and at large 
rapidities in H~b and J/y production covered by the LHCb experiment. 
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Figure 9: The differential cross sections do /dydpr of prompt J/\j/ production at LHC integrated over 
the specified y intervals. Solid, dashed and dashed-dotted curves correspond to the results obtained using 
the CCFM AO, BO and KMR uPDF. Dotted curves represent the contribution from sole direct production 
mechanism cal culated with the CCFM AO uPDF. The CMS data /^/. 
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Figure 10: The differential cross sections da j dydpj of the J/\j/ production at LHC integrated over the 



specified y intervals in comparison with the ATLAS data [34]. 
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Figure 11: The differential cross sections do / dydpr of the J/\j/ production at LHC integrated over the 
specified y inrervals in comparison with the LHCb data j^3$. 

Our study has demonstrated also that in the framework of the kj -factorization approach there is no 
need in a color octet contributions for the charmonium production at the LHC. 



As it was shown in [19] the future experimental analyses of quarkonium polarization at LHC are 



very important and informative for discriminating the different theoretical models. 
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